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Abstract 

Most discussions of chaotic scattering systems are devoted to two-dimensional systems. 
It is of considerable interest to extend these studies to the, in general, more realistic case 
of three dimensions. In this context, it is conceptually important to investigate the qual¬ 
ity of semiclassical methods as a function of the dimensionality. As a model system, we 
choose various three dimensional generalizations of the famous three disk problem which 
played a central role in the study of chaotic scattering in two dimensions. We present 
a quantum-mechanical treatment of the hyperbolic scattering of a point particle off a 
finite number of non-overlapping and non-touching hard spheres in three dimensions. We 
derive expressions for the scattering matrix S and its determinant. The determinant of 
S decomposes into two parts, the first one contains the product of the determinants of 
the individual one-sphere S-matrices and the second one is given by a ratio involving 
the determinants of a characteristic KKR-type matrix and its conjugate. We justify our 
approach by showing that all formal manipulations in these derivations are correct and 
that all the determinants involved which are of inhnite dimension exist. Moreover, for 
all complex wave numbers, we conjecture a direct link between the quantum-mechanical 
and semiclassical descriptions: The semiclassical limit of the cumulant expansion of the 
KKR-type matrix is given by the Gutzwiller-Voros zeta function plus diffractional cor¬ 
rections in the curvature expansion. This connection is direct since it is not based on any 
kind of subtraction scheme involving bounded reference systems. We present numerically 
computed resonances and compare them with the corresponding data for the similar two- 
dimensional 7V-disk systems and with semiclassical calculations. 
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1 Introduction 


Many if not most of the concepts in quantum chaos were developed and are being applied to 
two-dimensional systems. This is due to the relative simplicity of those systems as compared 
to three-dimensional ones. To find the periodic orbits in a generically three-dimensional ge¬ 
ometry, obviously requires much more work than performing the same type of analysis in 
a two-dimensional billiard, say. However, since the real world is three-dimensional, an ex¬ 
tension of chaos studies to more realistic systems is called for. The celebrated Hydrogen 
Atom in a strong magnetic field can be reduced to an, effectively, two dimensional system 
due to the cylindrical symmetry of the problem. Similar simplifications also exist in other 
systems. Recently, a full-fledged study of the three-dimensional Sinai-billiard has been per¬ 
formed [||. Here, we present a study of a chaotic scattering system in three dimensions. 
We investigate several generalizations of the two-dimensional three disk problem. We chose 
three-dimensional versions of this system since it played an important role in the development 
of many concepts and methods in chaotic scattering. In order to describe Wsphere scattering 
systems quantum-mechanically we extend the methods of Refs. [||, ^ to three dimensions. A 
related approach was already used in Ref. |Q] to describe the multiple scattering of an elec¬ 
tron in non-overlapping muffin-tin potentials. In Ref. |^] the scattering of a point particle on 
non-symmetric conhgurations of point-scatterers in three dimensions was investigated. 

In general, hyperbolic or even chaotic scattering systems have some advantages compared 
with bounded systems, if one is interested in the quantum-mechanical and the semiclassical 
description of a classically chaotic problem. Chaotic bounded systems are normally plagued 
by the presence of non-isolated, non-hyperbolic bouncing ball orbits (see, e.g., the Sinai- 
billiard or the stadium billiard), the importance of very long periodic orbits, and the problem 
that without fine-tuning a la Berry and Keating the semiclassical Gutzwiller-Voros zeta 
function |@, ^ would not predict real-valued energy eigenvalues. In contrast, the geometry 
of hyperbolic scattering systems can easily be chosen such that bouncing ball orbits are 
absent. Furthermore, the contributions of longer periodic orbits to the scattering matrix are 
automatically suppressed relatively to the shorter ones and no fine-tuning is necessary in 
order to predict scattering resonances as they are anyhow complex-valued. 

Two-dimensional iV-disk and three-dimensional iV-sphere systems are examples of such 
scattering systems which are simple enough to be studied in all detail semiclassically and 
quantum-mechanically. In the past the two-dimensional Sinai-billiard, which can be inter¬ 
preted as the scattering on an infinite, regular array of equal disks, has been quantized 
using the Korringa-Kohn-Rostoker (KKR) method |jl^. More recently the scattering of a 
point particle on hard disks in two dimensions |11] has been studied classically [12, 13], 
quantum-mechanically ||^ ^ and semiclassically 0, 0, |T|, 0 using the techniques of pe¬ 
riodic orbit theory |^. The range of validity of the purely geometrical semiclassical input 
has been investigated in Refs. |^, |^. In Refs.|0|, 0, ^-corrections to the geometri¬ 
cal periodic orbits were constructed, whereas the authors of Refs. 


24, 


extended the 

Gutzwiller-Voros zeta function to include diffractive creeping periodic orbits as well. Recently 
a formalism which also includes the limit of small disks (ka <C 1, k: wave number, a: radius of 
the disk) was presented 0] . In Ref. [^ the connection between the quantum-mechanical and 
semiclassical descriptions of V-disk scattering systems has been investigated in detail. On the 
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experimental side, the scattering on two equal disks was investigated using microwave cavi¬ 
ties [p^ . In Ref. the two-dimensional Sinai-billiard was treated in a scattering approach 
and in Ref. [29| diffractive effects were considered. 

Here we focus on the analogous scattering systems of N spheres in three dimensions. 
These systems are still simple enough to be treated quantum-mechanically, semiclassically 
and classically. The quantum-mechanical description of A^-sphere systems is similar to the 
two-dimensional A^-disk case, as essentially the same approach can be used. Also the semi- 
classical description with purely geometrical input is similar to the two-dimensional case. 
But diffractional corrections should be substantially different as the corresponding creeping 
orbits are now extrema on two-dimensional manifolds instead of one-dimensional ones. A de¬ 
tailed quantum-mechanical description of A^-sphere scattering systems is presented and the 
connection with the semiclassical treatment is investigated. 

The paper is organized as follows. In Sec. 2, we derive an explicit expression for the 
scattering matrix within the framework of stationary scattering theory using Green’s theorem. 
We do this in some detail to make the paper self-contained. In Sec. 3, the determinant of the 
scattering matrix is re-written as a product of an incoherent part and a coherent part. We 
show that the scattering resonances are given by the zeros of a KKR-type matrix. Moreover, 
we conjecture a direct link between the quantum mechanical and semiclassical descriptions of 
A^-sphere systems involving the Gutzwiller-Voros zeta function plus diffractional corrections. 
This connection is valid for all complex wave numbers. In Sec. 4, we give a proof that all formal 
manipulations performed in the preceding sections are allowed and that the determinants of 
the infinite matrices involved are well defined. In Sec. 5, we present numerical results on two-, 
three- and four-sphere scattering systems and compare them with the analogous two- and 
three-disk systems. In Sec. 6, we summarize our results and give an outlook. 


2 Calculation of the Scattering Matrix 

We describe the scattering of a point particle on N hard spheres within the framework of 
stationary scattering theory following the methods of Berry Q and Gaspard and Rice . In 
Sec. 2.1, we define the scattering matrix and outline our approach. In Sec. 2.2, the elements 
of the scattering matrix are worked out explicitly. 

2.1 Definitions and General Concepts 

To describe a generic configuration of N spheres we use the following notation: j £ {1,... , A^} 
specifies one spherical scatterer with radius aj. Rjji denotes the distance between the centers 
of the spheres j and j'. To specify positions we use A" -|- 1 coordinate systems: First we 
choose a global coordinate system (x, y, z) whose origin is situated at an arbitrary point in 
the neighbourhood of the N scatterers. This point is chosen to be the center of the large 
sphere of the integration volume used in Green’s theorem (see Sec. 2.2). If the scattering 
configuration possesses symmetries the origin of the global system is placed at the symmetry 
center of the entire system. In order to perform symmetry reductions we introduce N local 
coordinate systems, whose origins lie at the centers of the N spheres. The 

axes of these coordinate systems are chosen in such a way that the symmetry of the entire 
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configuration is respected. The vector from the origin of the global coordinate system to the 
center of the j-ih. sphere is called Sj and it is measured in the global system. All the vectors in 
the local coordinate system of sphere j are measured relative to this vector. /R~ji 

denotes the unit vector from the center of the sphere j to that of j' and it is measured in 
the (j)-system. In general, vectors with an upper index [j) are measured in the (j)-system, 
vectors without upper index are measured in the global system. 

A solution of the time-independent Schrodinger equation fulfills 

+ }?^ il){r) = 0, f outside the spheres, (1) 

'ip{f) = 0, r on the surfaces of the spheres. 

The energy of the particle is /2m and k is the wave vector of the incident wave. We 
expand the wave function ip{f) in a basis of eigenfunctions of angular momentum 

OO I 

W) = ^Lir)yi*m{k) ( 2 ) 

;=0 m=—l 


where k and k are the length and the solid angle of the wave vector, respectively. Because of 
this expansion, we construct solutions of the Schrodinger equation for the basis functions 

(V2 + p)V.L(r)=0. (3) 


Asymptotically for large distances from the scatterers {kr oo) the spherical components 
can be written as a superposition of in-coming and out-going spherical waves, 


*Lin - See 


oo V 


kr 


l'=0 m'=-l' 




(4) 


This equation defines the scattering matrix S. Its elements Sim,i'm' describes the scattering 
of an in-coming wave with angular momentum l,m into an out-going wave with angular 
momentum V. If there are no scatterers, we have S = 1 and the asymptotic expression of 
a plane wave is recovered. As is well known, the scattering matrix S is unitary because 
of probability conservation. 

To derive an explicit expression for the S-matrix, we use the free Green’s function, 

(V^ -|- P)G(r, f') = —47r(5^(f — r') , G{f,f') = G{f',f) (5) 

OO I 

G{f,f')=4mkY^ ji{kr<)hf\kr>)Yi*^{r)Yira{r) , ( 6 ) 

/=0 m=—l 
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where r< (r>) denotes the magnitude of the shorter (longer) of the two vectors r and r' 
whose directions can be expressed in any coordinate system. The Green’s function and the 
components wave function are inserted in Green’s theorem which yields 



'^L{r) (V2 + P)G(f, r') - G(f, r') (V^ + 

=-dvr / (fr'ijji^{r)5^{f-r') 

Jv 



da ■ 


^L(a)VG(a,f') - G{a,r'WUa) 
_ r 0 r' ^v, 

- 1-47rV>L(fO f'GG, 


(7) 


where V denotes the integration volume and dV is its boundary. The appropriate volume of 
integration V consists of a large sphere centered at an arbitrary point in the neighbourhood 
of the N scatterers. Its radius is chosen to be so large that asymptotic formulae like Eq. (Q) 
hold for points far away from the origin but still inside the integration volume. From this 
large sphere we exclude N spheres whose centers coincide with those of the N scatterers 
and whose radii are larger by an infinitesimal amount e than the corresponding radii of the 
scatterers. 


2.2 Computation of the Matrix Elements 

Equation (^), understood in the limit e —> 0, can now be applied to two different cases as we 
can choose the point f' to be outside or inside the volume V. In the first case f' is chosen to 
be some point on the surface of the scatterer j and therefore outside the integration volume 
V. In the second case r' is chosen to lie inside the volume V and its modulus r' is taken to 
be so large that asymptotic formulae like Eq. (Q) hold. The boundary dV decomposes into 
N+l disjoint parts: The outer boundary of the large sphere, dooV, and the N surfaces djV 
of the excluded spheres, which coincide with the scatterers in the limit e ^ 0. 

First case: f' = fj£ boundary of the j scatterer 
All components vanish on the surfaces of the scatterers, but the gradient of the wave 
function on these surfaces is nonzero. Its normal component can be expanded in spherical 
harmonics, defined on the surface of the scatterer, 

OO I' 

nj Y. ( 8 ) 

l'=0 m'=-l' 

where fij denotes a unit vector perpendicular to dV pointing outside V. The orientation 
is measured in the (j)-system. The unknown coefficients then uniquely 

determined through Green’s theorem (^). In this way the gradient of the wave function 
is characterized by the matrix which depends on the sphere label j and whose matrix 
elements are specified by angular momentum quantum numbers. Gollecting everything, we 
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arrive at a compact matrix formulation which expresses the gradient matrix A in terms of 
two matrices C and M 

a = A^' • . (9) 

where we use the Einstein-summation convention for the sphere indices j and jh The matrix 
elements of and M-^-^ read 


J oo I' 




and 


(47r) 

Zi=0m=-Z' 




X \J(21 + l)(2li + + 1) 


hi}\kaj) 




( 10 ) 


(h 

i' 


/ k 1' 


Vo 

0 

oj 

\m — m fa 

—m J 





’ oo V 

Zi=0 m=—l' 


-I 


X (21 + l)(2/i + + 1) 


( 11 ) 


[h 



( h 1' 


Vo 

0 

oj 

\m — fh fh 

—m J 


X Mka,) . 

h],’{kajf) 

Details are given in App. A. Here, ji{z) and h^^\z) are spherical Bessel functions and Hankel 
functions of first kind, respectively |^|. The vectors Sj are measured in the global coordinate 
system attached to the large sphere of integration and point from its origin to the center 
of the sphere (which is of radius Oj), Sj is its modulus, and Sj the corresponding unit 
vector. Rjjj is the vector from the center of the sphere j to the center of j' as measured 
in the local coordinate system attached to the j sphere, Rjji is its modulus and Rjji the 
corresponding unit vector. ^'^e the rotational matrices which 


and 

transform the local coordinate system of sphere j' to the global coordinate system and to the 


local coordinate system of sphere j, respectively. Finally, we use the definition of Ref. for 
the 3j-symbols. 

Second case: f'£V,r' large 

Because of Green’s theorem (|^ we have an explicit expression for V’z^(^) which yields inserted 
into the asymptotic expansion (0) an explicit formula for the scattering matrix 


S = 1 - i A^ • D-’ 


( 12 ) 
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where the elements of the matrix D-^ are found to be 

i 




jL„2 CXD 

S(-i) 


m+hjli-l 


( 13 ) 



h 


/ 1' h 

1 ] 

Vo 

0 

oj 

\m' rh — m' 

—rh 1 


X \J(21 + 1)(2V + l)(2/i + 1) 

Again, details are given in App. A. Here, the rotational matrix transforms the 

global coordinate system into the local coordinate system of the sphere. The gradient 
matrix A-^ appears again because of the boundary integrals on the scatterer surfaces. 

Hence, combining the results of both cases, we can use Eq. to eliminate the matrix 
from Eq. (IJ) |^]. We finally arrive at 


S = 1 - 


ia • ( m - 1 ) 


33 




(14) 


which is expressed in the global coordinate system. In order to shorten the notation we will 
suppress the labels j and /, unless otherwise specified. 

As mentioned above the effects of the scattering on the N spheres are described by the 
product —iC ■ M“^ • D which is of course the on-shell T-matrix. Examining the obtained 
formulae we see that M depends only on the relative positions of the scatterers and that it 
contains all the information about the geometry of the entire scattering system. It does not 
depend on the choice of the center of the large sphere in the integration volume and on the 
orientation of the corresponding global coordinate system. For these reasons we call M a 
characteristic matrix. It is of KKR-type |P, W|. In contrast, the matrices C and D depend 


on this choice and furthermore do not contain any information about the relative positions 
of the N scatterers. We may therefore conclude that the coherent multi-sphere part of the 
scattering is contained in M while the single sphere aspects are contained in C and D. 

If one would like to construct S explicitly from Eq. (14), it would be necessary to find the 
inverse of the infinite matrix M, which is a nontrivial task. But if one is only interested in 
spectral properties like scattering resonances it suffices to look for the poles of the determinant 
of the S-matrix |^, Q. The latter can be expressed in such a way that it does 

not involve any more (see the following Section). From Eq. ( p)4[ ) we expect that the 

resonances of the coherent part of the scattering are given by the zeros of the determinant of 

M. 


Before focussing on the determinant of S, we should have a closer look at Eq. (P) which 
is only determined up to a transformation of the form C’ = CE = AME = AM’. To get an 
unambiguous definition, we have chosen such a normalization that in the case of the scattering 
from only one sphere we have simply C = A. This choice implies M = 1 -|- W where the 
diagonal entries of W vanish. In Sec. 4 we sketch how to prove that Tr W converges absolutely 
in every basis so that the determinant of M is well defined [^]. This is very important as 
Det M plays a crucial role in calculating scattering resonances, as we will see now. 
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3 Determinant of the Scattering Matrix and the Link to the 
Semiclassical Zeta-Function 

Since we are interested in the scattering resonances it is sufficient to find the poles of the 
determinant, det S, of the scattering matrix [0, I 


as a function of complex wave 


number k. Obviously S is an infinite matrix, hence it is a nontrivial task to prove the existence 
of its determinant. To keep the discussion transparent, we postpone the formal proof to Sec. 4 
and anticipate the mathematical soundness of the calculations to be performed. We work 
out an explicit expression for det S in Sec. 3.1. In Sec. 3.2, we conjecture a direct link to the 
semiclassical zeta-function. 


3.1 Calculation of the Determinant 

A formal definition of the determinant of an infinite matrix Q = 1 + P is given by 

det(l + P) = exp{tr [ln(l + P)]} , (15) 

ln(l + P) = -Y.- -• (16) 

1 ^ 
n=l 

Equation (^), understood as Taylor-expanded expression (i.e., in the cumulant expansion, 
see below), is well defined if tr P converges absolutely in every basis |^. Formally one gets 
the following expression for the determinant of the A^-sphere S-matrix, 


det 


exp{tr ln(l — iC • M ^ • D)} 

{ CXD 

tr[(C-M-CD)"] 

n=l ^ 

{ jn 

- E —Tr [(M~i • D • C)*"] 

n ^ 

exp{Tr ln(l - • D • C)} 

Det (1 - iM-^ • D • C) 
Det[M-C (M-ID -C)] 

Det (M - iD • C) 

Det (M) ■ 


(17) 


Here, we introduced capital case traces and determinants Tr • • • and Det • • • in order to 
indicate that they refer to matrices labeled by the index triples j, I, m, whereas the lower case 
traces and determinants act on matrices which are just labeled by the (angular momentum) 
index pairs l,m. 






With the caveat that there might be also poles and zeros in the numerator, the resonances 
can be determined by just looking for the zeros of Det M in the complex /c-plane. In order to 
get a simpler expression for det which also gives more physical insight, we calculate the 
determinant of X = M — zD • C. With the knowledge of Eqs. (^) and (|l^ ) one can calculate 
the product D • C. In order to do this, the properties of consecutive rotations (changes of 
coordinate systems) and the re-coupling of angular momenta via 6-j-symbols (see Ref. |]^ 
for the necessary formulae) have to be considered. Here only the result will be given. 


D-? pJ 

,l'm' 


^ hY\kaj) 

oo V / \ ^ o 

q=0M=-i' V^i'/ * 

ji{kaj 


X ^{21 + 1){21' + l){2h + 1) IT"'' ji,ikRj,,)Yi,,m-Minify 

h), {kaji) 


V I \ f I' I 


0 0 0 / V M —m m — 


] ■ (18) 


It is easy to see that the matrix X = M — zD • C is given by 




i-iy 


j(i) 


'if) 


(19) 


where M* is a shorthand for (M(/c*))* and k) denotes the 1-scatterer S-matrix for the 

scattering from a single sphere with radius aj, described in a coordinate system whose origin 
lies in the center of the sphere, 


a' 


[kaji ) 

(kaji) 




( 20 ) 


The last equation is easily obtained by a comparison of the exact solution of this integrable 
problem and the ansatz (Q). The corresponding determinant, det A;), is of course 

independent of the coordinate system. We now have a (formal expression for the) product- 
structure for the determinant of X, 


DetX(A:) = Det(M(r))l det S(^)(j; fc) j 
Combining (|^) and (^), we finally obtain 


detS(^)(fe) 


Det (M(r))t 
Det M{k) 




( 21 ) 


( 22 ) 


The determinant splits up into an incoherent part, consisting of the product of the N de¬ 
terminants of the one-scatterer S-matrices, and a coherent part, given by the ratio of the 
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determinants of the hermitian conjugate of the characteristic matrix and the determinant of 
M itself. Equation (22) obviously respects the unitarity of the S-matrix, as the one-scatterer 


S-matrices ( |2^ ) are by themselves unitary because of {hm (ka))* = h!'^ {k*a) and as the co¬ 
herent part of detS^^^ is manifestly unitary. As mentioned, the scattering resonances are 
given by the zeros of Det M. However, this determinant does not only possess zeros but also 
poles. These poles cancel the resonance poles of the incoherent part of detS^^\ as DetM 
and the product of the det both involve the same number and power of Hankel functions 
h\^\kaj). The same is true for the poles of Det and the zeros of the product of the 
det both involve the same number and power of Hankel functions hf'\kaj). 

It is clear that det does not depend on the choice of the global coordinate system 
that was used in the definition of the S-matrix (Q) and to fix the center of the large spherical 
integration volume in Green’s theorem (j^. As we are interested in the coherent part of the 
scattering we will deal with DetM(A:) from now on. 

The symmetry of the scattering configuration leads to a block-diagonal form of the matrix 
M. If the symmetry group of the system is finite, we obtain [see App. B for the details of 
symmetry reductions] 


detS(^)(A:) 


]^detS(^)(j;A:)j 


OcCDet 


(23) 


where the index c runs over all conjugacy classes and Dc denotes the irreducible rep¬ 
resentation of dimension dc of the symmetry group. The last formula represents the final 
result of our formal treatment of scattering systems consisting of N spherical scatterers. The 
scattering resonances corresponding to the c-th irreducible representation of the symmetry 
group are given by the zeros of Det M 25 ^(A;) and they are dc-fold degenerate. 

In Sec. 4 it will be shown that all the formal manipulations that lead to Eqs. (p2[) and 


(^), and especially all the determinants appearing in it, are well defined if the number N of 
scatterers is finite and the scatterers do not overlap nor touch. 


3.2 Connection to the Semiclassical Gutwiller-Voros Zeta-Fnnction 

The determinant of the characteristic matrix M is understood in terms of the cumulant 
expansion |3|, which can be formally obtained from Eqs. (|I^) and (|l^), 


Det M 


exp{Tr logM} = exp{Tr log(l -|-W)} 


1 -t Tr W - - 
2 


Tr (W)2 - (TrW)^ 




(24) 


The first line of the upper equation is only of formal character — especially it is not defined 
at the zeros of Det M. Nevertheless, it offers an easy way to remember the second line, which 
defines the determinant of M, if W is trace class (see Ref. ||3| and references therein). 
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In a semiclassical description the scattering resonances can be extracted from the zeros 
of the Gutzwiller-Voros zeta function, which formally can be written as 01 


1 


ZGv{k-,z) = exp 




p r=i ^ I det(Jp'’ - i; 


(25) 


tp = exp 


u(m 


TT 

— Vt) — 
^2 


Here Jp denotes the monodromy matrix of the p-th primitive periodic orbit of topological 
length Up, Sp{k) and Vp are the corresponding classical action and Maslov index, respectively. 
The sum over r takes into account the repeated traversals of the primitive periodic orbits. In 
the above expression, z is a book-keeping variable and the Gutzwiller-Voros zeta function has 
to be evaluated aX z = 1. This expression is only of formal character, it has to be regulated. 


Expanding (^) in powers of z the curvature expansion is obtained, 


ZGv{k;z) = l-Z^6np,l 


f - 7 C 2 -h 

r p,r=l ^ 


C2 


p,r 


— 2 E E 


1 {hy 


p ^ 


r 






pLp' 


, 1 ,. . ,1 


p' 5-*- ■p ■p 

r p,r=ir p'^r=l 


(26) 

(27) 


= I det(Jp’' - 1)1 2 = |Ai,pA 2 ,pini - Ai,p-'’)(1 - A 2 ,p-n 


The n-th curvature contains all periodic orbits up to the topological length n. As can 


be seen from Eq. (^), the contribution of orbits of topological length n (> 1) gets reduced by 
pseudo-orbits, composed of shorter periodic orbits, of the same total length. Equations (p5|) 
and (^) (except for the last equality in the last line) are valid for two-dimensional V-disk 
systems as well as for three-dimensional V-sphere systems. In contrast to the well known 
two-dimensional case, in three dimensions the monodromy matrix, as being a 4x4-matrix, has 
two leading eigenvalues |Ai|, IA 2 I > 1 |^. For a generic periodic orbit in an V-sphere system 
the two leading eigenvalues are different. The most important exception is the only periodic 
orbit of the two-sphere-system, here Ai and A 2 coincide. In A^-sphere systems which have a 
two-dimensional analogue - this means that the centers of the N spheres are located all in one 
plane - one of the leading eigenvalues is given by the leading eigenvalue in the corresponding 
V-disk system while the other takes into account the instability of the periodic orbit against 
perturbations perpendicular to the plane. The lengths and Maslov indices of the periodic 
orbits in these V-sphere systems coincide with the corresponding quantities in the analogous 
V-disk system. As in the two-dimensional V-disk systems, orbits which are situated on the 


boundary of the fundamental domain require a special treatment [^, ^]. 
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It is easy to see that the series of Eq. (|24D has the same structure as the curvature ex¬ 
pansion of the semiclassical Gutzwiller-Voros zeta function (^). Furthermore, the quantum- 
mechanical description of A^-sphere scattering systems is analogous to the quantum-mechanical 
description of two-dimensional A^-disk systems. In the treatment of A^-disk systems a direct 
link between the quantum-mechanical and semiclassical descriptions has been established . 
On the quantum-mechanical side this link is based on {’M). Therefore we conjecture a similar 
direct link between quantum mechanics and semiclassics in A^-sphere systems: The semiclas¬ 
sical limit of the cumulant expansion in the Plemej-Smithies form (see Ref.j^]) of DetM is 
given by the curvature regulated Gutzwiller-Voros zeta function plus diffractional corrections, 


DetM(A:) 

Tr(W”(fc)) 


ZGv{k)\ curv. reg. 

^npr,n "-p 


(28) 


tpikY , 

rin-^ -h dittractional corrections. 


p,r 


Hence, the semiclassical limit of the n-th cumulant is given by the n-th curvature order in the 
Gutzwiller-Voros zeta function plus diffractional corrections. This direct link, which connects 
the quantum-mechanical and semiclassical descriptions, is valid for all complex k and not only 
for the isolated scattering resonance poles. We call it a direct link because it is not based on 
concepts characteristic to bound-state problems. So it is not based on any asymptotic limit 


of spectral densities (see Refs. []T| and 


hm (^N^^\k]b) - NY){k-b)) = ^ImTr lnS(A:) , 


(29) 


where (k; b) and (/c; b) are the integrated spectral densities belonging to two spherical 
bound systems, both of the same radius b, where one encircles the scattering region whereas 
the other doesn’t (see Ref. Q for a detailed discussion in the analogous two-dimensional N- 
disk systems). In fact, this formula is only correct in the double limit lim^^g where a 

small positive imaginary part ie has to be added to the wave number, before the limit 6 —> oo 
can be taken. 

In Ref. 139] a structure similar to that of Eq. (^) of the characteristic determinant for 
generic two-dimensional systems was derived in a semiclassical description under the Fred¬ 
holm theory. 


4 Justification of the Previous Calculations 

The derivation of the expression for the S-matrix (0) and the derivation of its determinant 
in Sec. 3 are of a purely formal character since all the matrices involved are of infinite size. In 
this Section, we show that the calculations of the previous Section are mathematically sound. 
This discussion is of considerable conceptual importance. However, those readers who are 
mainly interested in the numerical results of our study are adviced to skip the present Section. 

In the proofs that all the performed operations are well defined the so-called trace-class 
and Hilbert-Schmidt operators [^| play a central role. Trace class operators are those, 
in general, non-hermitian operators of a separable Hilbert-Space which have an absolutely 
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convergent trace in every orthonormal basis. An operator B belongs to the Hilbert-Schmidt 
class if is trace-class. Here we will not present all the proofs in detail. We refer 

to IQ] where the corresponding problem in the two-dimensional A^-disk scattering systems is 
treated. In that reference the most important properties of trace-class and Hilbert-Schmidt 
operators are listed: (i) any trace-class operator can be represented as the product of two 
Hilbert-Schmidt operators and any such product is trace-class; (ii) an operator B is already 
Hilbert-Schmidt, if the trace of B^B is absolutely convergent in just one orthonormal basis; 

(iii) the linear combination of a finite number of trace-class operators is again trace-class; 

(iv) the adjoint of a trace-class operator is again trace-class; (v) the product of two Hilbert- 
Schmidt operators or of a trace-class and a bounded operator is trace-class and commutes 
under the trace; (vi) if B is trace-class, the determinant det(l -|- zB) exists and is an entire 
function of z; (vii) the determinant is invariant under unitary transformations. 

In close analogy to the results of Ref. Qj the following steps can be proven - provided that 
N is finite and that the spheres do not touch nor overlap: 


(a) is a trace-class matrix for all complex k. Also is a trace-class matrix except at 

the isolated zeros of {kaj), where I is a nonnegative integer and j = 1,..., A^. (One 
very simple way to prove this is to transform D-^ and into the eigenbasis of the j 
sphere, see Eq.(pO|). In that eigenbasis both matrices become diagonal and the trace- 
class property can easily be checked by summing up the moduli of their eigenvalues. In 
the same way it can be checked that the one-sphere T-matrix k) — l is trace-class, 

too.) 

(b) Therefore the product D-^ C-^ is of trace-class as long as N is finite and is trace-class. 

(c) M-^-^ — 5^^'l = W-^-^ is trace-class, except at the same /c-values mentioned in (a). 

This can be proved by rewriting W-^-^ as the product of two matrices, and a 

diagonal matrix H-^-^ , which both can be shown — as in Ref.QI — to be Hilbert-Schmidt 
matrices: 


and 


G 


jj' 


(1 _ §n') ( ^ Jzi)!L 

V«T/ (2^ + 1)2 ^h!^J{kaaf) 

oo V _ 

X ^ ^ ihH'-y(2; + i)(2/i + i)(2/' + l) 

li=0 rh=—V 


(h 

i' 


/ h 1' 

^ ) 

Vo 

0 

oj 

\m — fh rh 

—m J 




(30) 


xrj 3 


= 5^ ^ 5vi"5„ 


{ 21 ' 


1)1 


3 \J (kctCLj 


h[}\kaji) 


(31) 
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where a > 2. This inequality is the reason that our proofs exclude the case of touching 
spheres. In fact, the geometry of the N scatters must be such that the inequality 
^aj + Qj' < Rjji is valid for all pairs j and j'. 

(d) Therefore M is bounded. 

(e) M is invertible everywhere where Det M is defined and nonzero (which excludes a 
countable number of isolated points in the lower /c-plane). Especially, M is invertible 
on the real k axis. Therefore, the matrix M~^ is bounded on the real k axis as well. 

(f) M-^DC, are all of trace-class (except at the isolated points mentioned 
in (a) and (e) as they are the product of a (finite number of) bounded and trace-class 
matrices, and tr [{0= Tr [(M-^DC) ”] exists (note the matrices on the 
l.h.s. are labeled by the index pairs l,m, whereas the ones on the r.h.s. are labeled by 
the index triples j, I, m). 

(g) M-iDC-1 is of trace-class because of (b) and (c) and the rule that the sum of two 
trace-class matrices is again trace-class. 


The properties (a) - (g) ensure that the derivation of Eq. (|l7|) starting from (H) is correct. 
Because of (e) formula ( p^ makes sense (also on the real A:-axis). Under the assumption that 
the matrix A, which characterizes the gradient of the wave function on the surfaces of the 
scatterers (see Eq. (|8|)), is bounded, the determinant of S on the basis of (p^) is defined. Now 
it is easy to see (using property (a) and (c)) that all determinants appearing in the formula 
(22) exists and that the unitary transformations leading to (23) are justified. 


5 Numerical Results 

We have calculated the scattering resonances of the three simplest and most symmetrical 
A-sphere systems. These systems consist of two, three and four hard spheres, respectively, 
which all have the same radius aj = a and the same center-to-center separation Rij = R: two 
non-touching spheres (two spheres), three spheres at the corners of an equilateral triangle 
(three spheres) and four spheres at the corners of a regular tetrahedron (four spheres). In all 
cases we have chosen the ratio R/a to have the fixed value 6 in order to be able to compare 
with older two and three-disk calculations. 

The quantum-mechanical resonances have been calculated as the zeros of the determinant 
of the characteristic matrix M in a finite basis corresponding to angular momenta from I = 0 
up to a certain maximum value which depends on the wave number k. We checked the 
accuracy of the results against a further enhancement of the basis. This method is applicable, 
as the trace class property of M — 1 guarantees the existence of the limit. The zeros in the 
complex fc-plane were determined with the help of a Newton-Raphson routine. 

The semiclassical resonances have been calculated as the zeros of the curvature regulated 
Gutzwiller-Voros zeta function without diffractional corrections. In the two-sphere and three- 
sphere systems the zeta function can be easily constructed under the techniques described in 
Sec. 3, with the input from the corresponding two-disk and three-disk systems, respectively. 
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The four-sphere-system is the simplest iV-sphere scattering system which does not have a two- 
dimensional analogue. In this case, only the three fundamental periodic orbits of topological 
length 1 have been determined so far. 

In Sec. 5.1, we discuss the quantum mechanically calculated resonances of the A^-sphere 
systems. They are compared to those of the two-dimensional A^-disk systems in Sec. 5.2. In 
Sec. 5.3, we compare the quantum mechanical and the semiclassical results for the A^-sphere 
systems. 

5.1 Quantum Mechanically Calculated Resonances 

In Secs. 5.1.1, 5.1.2 and 5.1.3 we discuss the two-, three and four-sphere systems, respectively. 

5.1.1 Two Spheres 

The scatterer configuration is shown in Fig. 1. This system has a continuous symmetry, 
the rotational symmetry about the axis which joins the centers of the spheres. Hence, the 
corresponding symmetry group, D^ch, is infinite and has 4 one-dimensional irreducible repre¬ 
sentations and an infinite number of two-dimensional ones |^. From the 4 one-dimensional 
representations only two are present due to the lack of any inner structure of the two spherical 
scatterers. 

The blocks of the symmetry reduced M-matrix are given by: 
one-dimensional representations (with m = 0): 

M„.(D) = sw + + mi' + 1) 

hj^, (ka) 

+ !)(('' g)Yhj;>(kn) 

li=0 

^ r 0 for Dc = Aig 
1 1 for Dc = Aiu 

two-dimensional representations: 

Mh'(^) = + (~1)'^ (~1)* (2^ + 1)(2/' -|- 1) 

hi, {ka) 

^ )h\l\kR) 

m —m y " 

^ r 0 for He = Hmg 
\ 1 for Dc = Emu 

1,1' > m > 0 , m fixed. 
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The degeneracy-degree of the resonances coincides with the dimensionality of the repre¬ 
sentation. As a consequence of the continuous rotational symmetry, \m\ is a good quantum 
number and there are two irreducible representations corresponding to each value of \m\. 

The numerically calculated resonances of the two-sphere-system are shown in Figs. 4-6. 
The leading resonances of Fig. 4 form a regular structure which we call a Gutzwiller-band. 
The suppressed resonances do not build up such a regular structure; we say they lie in 
a diffraction-band. This terminology is chosen in analogy to the similar results found in 
two-dimensional problems (scattering from N hard disks) I, I, [^ , [^ . The resonances 
corresponding to increasing \m\ become increasingly suppressed. This is due to the centrifugal 
barrier, which in cylindrical coordinates reads vr? jr^. 

In Fig. 5 we see that the leading resonances of the \m\ = 1 representations for small real 
parts of k lie in the diffraction band. For real parts of k bigger than approximately 10/a two 
Gutzwiller bands build up, which are shifted by half a spacing in Re k but coincide in their 
values of Im k. 


5.1.2 Three Spheres 

The scatterer configuration of the three-sphere-system is shown in fig 2. The symmetry 
group D'ih consists of 12 elements. It has four one-dimensional and two two-dimensional 


representations |^]. 

The symmetry reduced expressions for the distinct blocks of the M-matrix are given as 
follows. 

There are four one-dimensional representations (0 < m, m')'. 


(^1 ) 

M/m,/'m'(^2^0 




(a = 0) 


0 

E 

0 

~ 1 ) 

0 

0 


E 




{a = 0) 


{I + m) or {V + m') odd 
otherwise 

{I + m) or {V + m') even 
otherwise 

{I + m) or {V -|- m!) odd 
otherwise 

{I + m) or {V + m!) even 
otherwise 


(32) 

(33) 

(34) 

(35) 
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In addition, there are two two-dimensional representations (m, m' integer): 

(36) 

(37) 


with 

h)}\ka) 

X f;i'V(2/ + l)(2r + l)(2/i + l)(jj Q L 

/i=0 

X h\l\kR)Yi^^rn.-m'(^,0) cos ^^(5m - m'- 4)^ . 





0 

0 


{I + m) or {I' + m!) odd 
otherwise 

{I + m) or [I' + m!) even 
otherwise 


The numerically calculated resonances of the Ai and Ai" symmetry classes are shown in 
Fig. 7. The Ai" resonances are suppressed compared with the Ai resonances, because the 
wave function has to vanish in the plane of the triangle due to the symmetry (see character 
table of D^fi in Ref. ^^). This is the plane which contains all geometric periodic orbits 
of the three-sphere-system, which play the dominant role in a semiclassical treatment in the 
considered wave number regime. The resonances of the A 2 \ A^' and E', E" symmetry classes 
show a similar behaviour. 
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5.1.3 Four Spheres 

The scatterer configuration of the four-sphere-system is shown in fig 3. Compared with the 
two- and three-sphere-systems this is the first genuine three-dimensional scattering system 
since the scatterers (and therefore the set of all geometric periodic orbits) do not lie in a 
plane. 

The symmetry group has 24 elements. It has 2 one-dimensional, 1 two-dimensional 
and 2 three-dimensional irreducible representations 1^^. Here only the one-dimensional rep¬ 
resentations are considered. The corresponding blocks of the symmetry reduced M-matrix 
are given by: 




Q 


9m 


dlnm'iP) 

cos(6»o) 

cos(/Io) 




(38) 


-m-l ji{ka) 


T „ \/47ri 9m9m' 


2 ■ ■■■' 


h\, {ka) 


OO I' 


xE E i^hf\kR)^{2l + l){2l' + l)(2r + 1) 

[=0 M=-l' 

X ( — 1 )^ (dln'MiPo) + ( — 


I V I 
0 0 0 


(-i)”n,„-«(«o.o)( 

+ (-i)‘’n,-„-M(»o.o)( 


I I' I 

m — M M —m 


I V I 

—m — MM m 


0 for Z) = Ai 
1 foT D = A 2 

\/2 for m = 0 
1 for m = 3, 6,9,..., Z 
0 otherwise 

{jm 1 \jm') 


2 

To 

1 

3 


sin(6'o) = 


1 


v/3 


sin(/Io) = -V 2 


The numerically calculated Hi-resonances are shown in Fig. 8. 

A comparison of the first leading Gutzwiller resonances of the two-, three- and four- 
sphere-systems shows that the spacing in the real part of k in the resonance band slightly 
decreases if one looks hrst at the two-sphere, then at the three-sphere and finally at the 
four-sphere system. This spacing is governed by the inverse length of the (averaged) periodic 
orbits of topological length one in the fundamental domain. In the two-sphere case there 
exists only one periodic orbit. The three-sphere system has in addition to this orbit one 


18 






further fundamental periodic orbit, which in the global domain corresponds to an equilateral 
triangle spanned inbetween the three scatterers. The length in the fundamental domain of 
this second periodic orbit is only slightly bigger than that of the two-sphere orbit. In the 
four-sphere system there are three fundamental periodic orbits: The two orbits of the three- 
sphere system and an additional orbit which touches all four spheres in the global domain and 
which again it slightly bigger than the other two when measured in the fundamental domain. 
Thus the average length of the fundamental orbits increases with the increasing number of 
spheres. 

As far as the imaginary parts of the first leading Gutzwiller resonances of the two-, three- 
and four-sphere-systems are concerned, we see that the two-sphere resonances are slightly 
more suppressed than the three-sphere resonances, which are in turn slightly suppressed as 
compared to the four-sphere-resonances. This is due to the fact that the addition of one 
further sphere increases the probability that the particle is rescattered to the other spheres 
and therefore trapped for a longer time in the scattering region. 

5.2 Comparison of the Quantum Mechanically Calculated Resonances of 
Al-Sphere and Al-Disk Systems 

If the centers of the spheres of a given iV-sphere configuration lie all in one plane, there 
exists an analogous two-dimensional iV-disk scattering system. In this case it makes sense to 
compare directly these analogue systems, although they differ in their dimensionality: they 
have in common the entire set of classically allowed periodic orbits which play the dominant 
role in a semiclassical description in the wave number regime considered here ]l^ , |^ . 

In Fig. 9 the resonances of the totally symmetric representations of the two-sphere and 
two-disk systems are shown. In both cases the leading resonances lie in Gutzwiller bands. 
Suppressed resonances form diffraction bands. However, in the two-disk system a sub-leading 
Gutzwiller band seems to appear for larger wave numbers (Re A: > 15/a). The leading 
resonances in the two-sphere and two-disk systems have the same real part, but the whole 
two-sphere band is shifted down into the negative complex k plane, because the only existing 
geometrical periodic orbit is more instable in the three-dimensional case than in the two- 
dimensional one. In the semiclassical formulae the higher instability is easily explained by 
the existence of a second leading eigenvalue of the monodromy matrix (see Sec. 3.2). There 
is no such coincidence in the real parts of the resonances of the diffraction bands. This 
behaviour is also expected from a semiclassical point of view as the diffractional orbits in 
both cases act on different manifolds. 

In Fig. 10 the resonances of the totally symmetric representations of the three-sphere 
and three-disk systems are shown. In both cases the leading resonances lie in Gutzwiller 
bands. Suppressed resonances form diffraction bands. As in the two-scatterer systems, we 
find a good agreement between the real parts of the sphere and disk resonances (which is 
getting better for larger wave numbers) while the whole three-sphere Gutzwiller band is 
shifted to smaller imaginary parts by the same amount as in the two-scatterer case. The 
same behaviour is observed for the two-dimensional representations E' (three-sphere) and E 
(three-disk), respectively. 

In summary we see, without using the results of any semiclassical calculation, that the 
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leading resonances in the analogous two- and three-dimensional A^-scatterer systems are dom¬ 
inated by the contribution of the geometric periodic orbits. In contrast, the suppressed res¬ 
onances show no one-to-one correspondence between the disk and sphere resonances. So the 
conclusion, which so far is only based on the presented quantum-mechanical data, is that the 
suppressed resonances are due to diffractional effects, which should be different in two and 
three dimensions. 


5.3 Comparison of Quantum Mechanically and Semiclassically Calculated 
Resonances in A^-Sphere Systems 

As described in Sec. 3.2, it is possible to transfer the well developed techniques of calculating 
scattering resonances under the periodic orbit approximation with only classically allowed 
orbits 1111,111 from A^-disk systems to scattering problems with N hard spheres in three 
dimensions. 

In particular, the semiclassical resonances of the two-sphere-system are given by the zeros 
of the following spectral zeta function. 


g/u 


z{k) = n 

c=0 


\m\ 


{k) 


00 / ^ik{R-2a) \ 


(39) 


where the ± sign refers to the g/u representations and contains all terms of a given 
\m\ in the last expression. Aq is the leading eigenvalue of the monodromy matrices of the 
2-disk and the two-sphere systems. In the three-dimensional case this eigenvalue is two-fold 
degenerate. 

Also in the three-sphere system the semiclassical zeta function can be constructed from 
the corresponding expression for the two-dimensional three-disk system. As all periodic orbits 
lie in one plane, the 4x4 monodromy matrix decomposes into blocks of 2 x 2 matrices, where 
the off-diagonal blocks vanish. One block describes the motion in the plane and it is given 
by the monodromy matrix of the three-disk system. The other block describes the motion 
perpendicular to the plane and can be constructed from the same section lengths k, local 
curvature radii pi and scattering angles Oi as in the three-disk system. The only difference is 
that the matrix elements of the 2x2 dimensional reflecting matrix of Ref. [IT^] change to |R| 


/ -I -2/[piCos(6»i)] 

I 0 -1 


/ 1 2cos(6'i)//9i \ 

lo 1 ; 


whereas the translational matrices stay as 



(40) 


Note that in the 2-scatterer cases the angles Oi are always zero. 

In Fig. 11 the quantum-mechanically and semiclassically calculated resonances of the com¬ 
pletely symmetric representation of the two-sphere system are shown. The agreement between 
the quantum-mechanical and semiclassical data in the leading Gutzwiller band is very good 


20 




except for the first few resonances located at small real parts of the wave number k. In 
contrast, the suppressed quantum resonances situated in the diffraction band cannot be de¬ 
scribed by these semiclassical calculations. Note that the agreement of the leading resonances 
calculated quantum-mechanically and semiclassically is already very good at quite small real 
parts of k. For increasing Re A: this agreement improves. 

Also the quantum-mechanical and semiclassical resonances of the completely symmet¬ 
ric representation of the three-sphere-system have been calculated. Again, we find a good 
agreement between the quantum-mechanical and semiclassical data in the case of the leading 
Gutzwiller band except for the first few low-lying resonances. As in the two-sphere case 
the suppressed quantum resonances of the diffraction band cannot be described by these 
semiclassical calculations. 

In the four-sphere-system, we have only determined the three fundamental periodic orbits 
so far. Even in this case a comparison between the zeros of the quantum determinant Det M 
calculated via the cumulant expansion (^) which has been truncated after the first cumulant 
shows an agreement between the quantum and semiclassical leading resonances. 

6 Conclusions 

We presented a quantum mechanical and a semiclassical discussion of a chaotic system in 
three dimensions. As a model, we chose the scattering of a point particle on N hard spheres 
in three dimensions. We showed that certain methods which were developed in the frame¬ 
work of two-dimensional systems can be extended to three dimensions in a straightforward 
way. Within the framework of stationary scattering theory and using Green’s theorem an 
expression for the determinant of the scattering matrix S was derived. In order to determine 
the scattering resonances it suffices to look for the poles of the determinant of the S-matrix. 
The determinant of the S-matrix of the entire Wsphere scattering system splits up into a 
coherent and an incoherent part. The coherent part is given by the ratio of the determinants 
of the hermitian conjugate and of the characteristic matrix M itself, where M contains the 
information on the geometry of the entire A^-scatterer configuration. The incoherent part 
consists of the product of the N determinants of the 1-scatterer S-matrices, each expressed 
in a coordinate system whose origin lies in the center of the individual scatterer. The ex¬ 
pression for det respects the unitarity of the scattering matrix and guarantees a unitary 
semiclassical limit without any resummation techniques a la Berry and Keating |^]. The 
scattering resonances are given by the zeros of the determinant of M. The poles of the inco¬ 
herent part of det get canceled by the poles of Det M. A proof was given that all the 
formal manipulations in our derivations are allowed and that the final expression for det 
is well defined provided that the number of scatterers is finite and that they do not touch 
nor overlap. Similar results have been found in the description of two-dimensional Wdisk 
scattering systems [^, 

We conjecture the following direct link between the quantum-mechanical and semiclassical 
descriptions of Wsphere scattering systems: The semiclassical limit of Tr (W”), W = M — 1, 
is given by those terms of the curvature regulated Gutzwiller-Voros zeta function which 
correspond to periodic orbits of total topological length n = UpV (r denotes the number of 
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repeats), each weighted with the topological length Up of the underlying primitive periodic 
orbit, plus diffractional corrections. As the determinant of M is given by the cumulant 
expansion, this means that the semiclassical limit of the n-th cumulant of Det M is given by 
the n-th curvature order of the Gutzwiller-Voros zeta function plus diffractional corrections. 
This connection holds for all complex wave numbers. It is direct, as it does not rely on the 
subtraction of the (integrated) spectral densities of two equally sized infinitely large bounded 
reference systems - one containing the A^-sphere scatterer and the other not 
this is analogous to the results found in A'-disk scattering systems [^. 

Qualitatively, the distribution of the quantum-mechanically calculated resonances of N- 
sphere systems is similar to the known results of A-disk systems § i- Comparing the 
computed resonances of analogous A-sphere and A-disk systems, which have all the geomet¬ 
ric periodic orbits in common, we see that the leading Gutzwiller resonances have the same 
real part but the whole resonance band of the three-dimensional system is shifted to smaller 
imaginary parts. This is due to the increased instability of the geometrical periodic orbits 
in the three-dimensional case. In the case of the suppressed resonances located in diffraction 
bands there is no such correspondence. The comparison of quantum-mechanically and semi- 
classically calculated A-sphere resonances shows qualitatively the same behaviour as in the 
two-dimensional A-disk systems. 

The semiclassical investigation of A-sphere scattering systems is still a sparsely studied 
field. The description of diffractional corrections could be the subject of future projects. In 
this context an investigation of A-sphere systems in which the sphere separation is much 
bigger than the radii of the scatterers would be of interest, too, because diffraction effects 
should become more important. Also the opposite case with configurations of almost touching 
scatterers provides an interesting system in which bound states might develop. Another 
important point is a further investigation of the here conjectured link between the quantum- 
mechanical and semiclassical descriptions. 


17 , 


Again 


Acknowledgment 

We would like to thank Friedrich Beck and Predrag Cvitanovic for fruitful discussions and 
helpful advice. We are grateful to Harel Primack for sending us the structure of the three- 
dimensional A-sphere Jacobians. A.W. acknowledges the warm hospitality of the Genter of 
Chaos and Turbulence Studies at the Niels Bohr Institute in Copenhagen, Denmark, where 
part of this work was done. T.G. thanks the Danish Research Council for financial support. 

A Details of the Calculation of the Scattering Matrix 

We distinguish between the two possible cases defined in Sec. 2.2. 
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First case: r' = rj £ boundary of the j scatterer 
Green’s theorem can be written in the form 


N 

Q = iL + T. If ^ 

j'=i 


( 41 ) 


where the integrals are given by 


P = 


WooV 


da - [V’L(«)VG(a,F,) - G{a,fj)V'tpi^{d)] 


4 = - 


'dj,V 


da ■ G{a,fj)V'tl}f^{a) . 


(42) 

(43) 


The integrals defined in Eqs. (42) can be worked out in a straightforward calculation, we 
arrive at 


li, = ^ (-l)™z^i+^y(2Z + l)(2Zi + l)(2/2 + l) 

h,l2=0 fh,m2=-l2 


X JhikspjiPkaj) 


h h I 


h 


h I 


0 0 0/ \m — 777-2 7772 — 777 

X yh,m-m2{Sj)Yi2,fh{df)D^l^2^gl,j) , 

oo I' 

Jj = 4TTikJ2 H (^‘j^irn,l'rn'P(kaj)h\}\kaj)Yi,rn'{a^j^) , 

l'=0 m'=-l' 


(44) 

(45) 


p = 

7' 


oo I' h 

^^ka} E E E 

^^^l,^ 2=0 m'=—l' m 2 ,ih=—l 2 


h h I' 


h h I' 


0 0 0 ) \m' — 7772 7772 — 777^ 


X \J(2U + l)(2/i + 1)(2Z2 + 1) 

X p {ka,,)ji, {kap {kR,,pYi,^^,.^, {R^]) Yi,^{df) (j', j) . (46) 
In the derivation we made use of the addition theorems for Bessel functions, 


OO ll 

i’'jl{kr)Yim{r) = E E (“I)”"* 

ll yl2—0 mi=—li 




X \J (2Z + l)(2Zi + 1 )( 2/2 + PJh (ks) ji^ (ka) 

), (47) 


h h l\( h h I 

0 0 0 / \mi m — mi —m 


r = s + a , 


23 






and Hankel functions 


i^h\^\kr)Yim{f) 


r = 


h,h=0 mi=—li 

X (2/ + l)(2/i + 1)(2/2 + {ks)ji^{ka) 

f n n ^ ) y/,™,(s)ll,,„,(a), (48) 

\ 0 0 0 / \mi m — mi —mj i- i-'- ' 

s + a , s > a , 


which can be easily proven using = e**^'*e**^'“ and the known properties of Bessel and 


Hankel functions |^. We use the definition of Ref. [^] for the 3j-symbols. 

For large distances from the scatterer conhguration we use the following asymptotic ex¬ 
pressions for Hankel-functions (kr 


oo 




h\ >{kr) 


kr 


h\^\kr) 


J_g+i(A:r-i±i7r) 

kr 


(49) 


In changing coordinate systems we adopt the definitions of Rose 
and the irreducible representations of the rotational group, e.g., 


for the Euler angles 




and Dl^,^{gl, j) denotes the corresponding quantity where the Euler angles describe the ro¬ 
tation of the axes of the (j)-system in those of the global system. For the spherical harmonics 
Yim we adopt the usual definition. 

In particular, in the calculation of the integration is performed in the global coordinate 
system. The Green’s function are evaluated under the asymptotic expression ( |4^ ) and finally 
the addition theorem (^7|) is applied and the coordinate system is changed. This leads to 
Eq. (^^. In the calculation of /j the integration is performed in the (j)-system. One easily 
obtains Eq. (45). In the calculation of Ij,, j ^ f, we hrst perform the angular integration in 
the (j^)-system. Then the addition theorem (|4^ is used and finally the coordinate system is 
changed. We get Eq. (^). 

The expressions {M) to (^) are already written in a form R is 


I), written as represents an equality between 


now easy to see that Eq. 
functions defined on the surfaces of the j-th scatterer and that it becomes 


E CL,l,m,yhm,{af) = Y 








(50) 


l2,m2 


j' yl' ,Tn' 


We have obtained an equality between the coefficients appearing in the last equation up 
to a transformation of the form C’ = CE = AME = AM’. Because of the reasons already 
discussed in Secs. 2 and 4 we write the characteristic matrix as M = 1-|-W where the diagonal 
entries of W vanish. With these definitions we obtain Eqs. §)-&■ 
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Second case: f'£V,r' large 
According to we obtain 

&(-') =-A f/i' + f/fj 

where the integrals are given by 


(51) 


r = 

oo 


Id^v 


, da • (^L(a)VG(a, f') - G(a, r' )WL(5)) 


(52) 

(53) 


which yields Eq. (^) and thus determines the components of the wave function. The evalu¬ 
ation of the integrals yields 


/; = - / da-G{a,r'WUa) 

JdjV 


Ilo = ji{kr')Yim{r) 


(54) 


r 

j 


oo I' oo h 

= {4TT)hka]Y. Y. (-1) 

m'=—l' li ,/ 2=0 mi=—li 


Tfl' -\-l2 ^^1+^2 — 


X {21' + l){2li + l){2l2 + l)hf^ {kr)ji^ {ksj)ju (kaj) 


h h I' 
0 0 0 


h 


/ / ) A/o.m'—mi (Si)!^! mi (f* ')-^lml'm' ' 

mi m — mi —m ' ’ i iv / iiu,i m 


(55) 


The calculation of can be carried out in an analogous way as the determination of 
In the calculation of /J the integration is performed in the (j)-system, then the addition 
theorem (|4^) is used. This leads to Eq. (^). 

One could insert the last results in Eq. (51) and drop the condition kr' ^ 1. This would 
yield an expression for the wave function valid for all r'. The only restriction to r' is that the 
addition theorem (|4^ of Hankel-functions has to hold. In order to determine the S-matrix 
it suffices to use r' far away from the scatterers, such that Eq. (^) becomes valid. Hence, 
inserting (^) in (^) we obtain 

^lm,l'm' — dllidn 


- E 


■*-*Z"m" ,l'm' 


(56) 


We use A' 


Iml'm' 


— spl 
L^ih=—l 


^im 9^) result Compatible with the results 


obtained in the first case. By this operation the N local coordinate systems appear and we 
obtain Eqs. (^) and ([T^ . 


B Symmetry Considerations 

In order to simplify the determination of scattering resonances, given as the zeros of the 
determinant of the characteristic matrix M, the symmetry of the scattering configuration 
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can be used to block-diagonalize M. (These operations are allowed as M — 1 has been proven 
trace-class and as Det M exists.) A similar approach was used in Ref. [^] in order to perform 
symmetry reductions in A^-disk scattering systems in a semiclassical description. 

The infinite square matrix M can be interpreted as a linear application acting on the 
square integrable functions defined on the surfaces of the N scatterers. In Sec. 2 we have seen 
that M contains the complete information on the geometry of the entire scattering system, 
so we have 

[M,p(<7)]=0, ygeG. (57) 

Here p denotes the linear representation of the symmetry group G acting in the same space 
as M. As the entire scattering configuration can be constructed by successive applications 
of the symmetry transformations of G acting on an appropriate fundamental domain, the 
representation p is decomposable, 

p = R® p . (58) 


The representation R acts on the IGI (= number of elements of G) copies of the funda¬ 
mental domain and p acts on an appropriate basis of the square integrable functions defined 
in the fundamental domain. Now it is easy to see that R is the regular representation |40| 
of G. If the symmetry group G is finite, the regular representation is, as is well known, 
reducible, 

R = dcDc ■ (59) 


Here the index c runs over all conjugacy classes and Dc is the c-th irreducible representation of 
dimension dc- Therefore, if we are dealing with scattering systems that possess only discrete 
symmetries (this suffices to ensure |G| < oo as we are dealing with finite, N < oo, systems) 
we obtain 

Det M = ]^(Det ■ (60) 


The last result is inserted in Eq. (|^) and formula (^) is obtained. 


26 



References 


[ 1 ] 

[ 2 ] 

[3] 


H. Primack and U. Smilansky, Phys. Rev. Lett. 74 (1995), 4831-4834. 
P. Gaspard and S. a. Rice, J. Chem. Phys. 90 (1989), 2255-2262. 


A. WiRZBA AND M. Henseler, preprint, September 1996, 


see: http://crunch, ikp.physik.th-darmstadt.de/~wirzba/ or 
[http : //www. nbi. dk/^predrag/QCcourse/|. 


[4] P. Lloyd and P. V. Smith, Adv. Phys. 21 (1972), 69-142. 

[5] W. John, B. Milek, H. Schanz and P. Seba, Phys. Rev. Lett. 67 (1991), 1949-1952. 

[6] M. V. Berry and J. P. Keating, Proc. Roy. Soc. A 437 (1992), 151-173. 

[7] M. C. Gutzwiller, “Chaos in classical and quantum mechanics,” Springer (New York), 
1990. 


[8] A. VOROS, J. Phys. A 21 (1988), 685-692. 

[9] M. V. Berry, Ann. Phys. (N.Y.) 131 (1981), 163-216. 

[10] J. Korringa, Physica 13 (1947), 392-400; W. Kohn and N. Rostoker, Phys. Rev. 
94 (1954), 1111-1120. 

[11] B. Egkhardt, .j. Phys. A 20 (1987), 5971-5979. 

[12] P. Gaspard and S. A. Rice, J. Chem. Phys. 90 (1989), 2225-2241. 

[13] P. Gvitanovic, “The power of chaos,” in “Applications of Chaos, Electrical Power 
Research Institute workshop 1990” (J. H. Kim and J. Springer Eds.) John Wiley &: 
Sons, New York 1992. 

[14] P. Gaspard and S. A. Rice, J. Chem. Phys. 90 (1989), 2242-2254. 

[15] P. Gvitanovic and B. Eckhardt, Phys. Rev. Lett. 63 (1989), 823-826. 

[16] P. Gvitanovic, P. E. Rosenqvist, G. Vattay and H. H. Rugh, CHAOS ^ (1993), 
619-636. 

[17] B. Egkhardt, G. Russberg, P. Gvitanovic, P. E. Rosenqvist and P. Scherer, 
in “Quantum chaos between order and disorder” (G. Casati and B. Chirikov Eds.), 
Cambridge University Press, 1995, pp. 405-433. 

[18] A. WiRZBA, Nucl. Phys. A560 (1993), 136-150; CHAOS 2 (1992), 77-83. 

[19] P. Gvitanovic, G. Vattay and A. Wirzba, “Quantum Fluids and Classical Deter¬ 
minants,” to appear in ’’Classical, Semiclassical and Quantum Dynamics in Atoms - 
in Memory of Dieter Wintgen,” (H. Friedrich and B. Eckhardt Eds.), Lecture Notes in 
Physics (Springer Verlag, Berlin). 


27 







[20] P. Gaspard and D. Alonso, Phys. Rev. A 47 (1993), R3468-3471; 
D. Alonso and P. Gaspard, CHAOS Z (1993), 601-612. 


[21] P. Gaspard, in “Quantum chaos between order and disorder,” (G. Casati and B. 
Chirikov Eds.), Cambridge University press, 1995, pp. 385-404. 


[22] G. Vattay, |chao-dyn/9406005| (1994); 

G. Vattay and P. E. Rosenqvist, Phys. Rev. Lett. 76 (1996), 335-339. 


[23] G. Vattay, A. Wirzba and P. E. Rosenqvist, Phys. Rev. Lett. 73 (1994), 2304- 
2307. 


[24] G. Vattay, A. Wirzba and P. E. Rosenqvist, Proceedings of the International 
Conference on Dynamical Systems and Chaos, Tokyo, May 1994 (Y. Aizawa, S. Saito 
and K. Shiraiwa Eds.), World Scientihc, Singapore 1995, Vol. 2, pp. 463-466. 

[25] P. E. Rosenqvist, G. Vattay and A. Wirzba, J. Stat. Phys. 83 (1996), 243-257. 

[26] P. E. Rosenqvist, N. D. Whelan and A. Wirzba, J. Phys. A 29 (1996), 5441-5453. 

[27] A. Kudrolli and S. Sridhar, Proceedings of the 4th Drexel conference (World Scien¬ 
tific, to appear), see: ^ittp : //sagar.physics.neu. edu/preprints. htmlj (1995). 

[28] H. Schanz and U. Smilansky, Chaos, Solitons and Fractals 5 (1995), 1289-1309. 

[29] H. Primack, H. Schanz, U. Smilansky and I. Ussishkin Phys. Rev. Lett. 76 (1996), 
1615-1618. 


[30] M. Abramovitz and I. A. Stegun (Ed), “ Handbook of Mathematical Functions with 
Formulas Graphs and Mathematical Tables,” Dover Pub. Inc. New York, 1964. 

[31] Landolt-Bornstein, “Zahlenwerte und Funktionen aus Naturwissen- 

schaften und Technik, Band 3: numerical tables for angular correlation computations in 
a-, (3-, and 7 -spectroscopy: 3j-, 6 j-, 9j-symbols, F- and P-coefficients,” Springer-Verlag, 
1968. 

[32] P. Lloyd, Proc. Phys. Soc. 90 (1967), 207-216. 

[33] M. G. Krein, Mat. Shorn. (N.S.) 33 (1953), 597-626; Sov. Math.-Dokl. 3 (1962), 707- 
710. 

M. S. Birman and M. G. Krein, Sov. Math.-Dokl. 3 (1962), 740-744. 

[34] J. Friedel, Nuovo Cim. Suppl. 7 (1958), 287-301. 

[35] M. Reed and B. Simon, “Methods of Modern Mathematical Physics, Vol. IV: Analysis 
of Operators, Chap.XIII.17,” Academic Press (New York), 1976; 

B. Simon, Adv. Math. 24 (1977), 244-273. 

[36] B. Lauritzen, Phys. Rev. A 43 (1991), 603-606; 

J. M. Robbins, Phys. Rev. A 40 (1989), 2128-2136 


28 







[37] P. CviTANOVic AND B. ECKHARDT, Nonlinearity 6 (1993), 277-311. 

[38] R. Balian and C. Bloch, Ann. Phys. (N. Y.) 85 (1974), 514-545. 

[39] B. Georgeot and R. E. France, Phys. Rev. Lett. 74 (1995), 2851-2854, 4110-4113; 
S. Fishman, B. Georgeot and R. E. France, J. Phys. A 29 (1996), 919-937. 

[40] M. Hammermesh, “Group Theory and its applications to physical problems,” Addison- 
Wesley, 1962. 

[41] H. Primack, private communieation. 

[42] M. E. Rose, “Elementary Theory of Angular Momentum,” John Wiley and Sons Inc. 
New York, 1957. 


29 




Figure 1: The two-sphere-configuration consisting of two equal spheres. The global and local 
coordinate systems used are also shown. 



Figure 2: Quantum-mechanically calculated Ai^-resonances of the two-sphere system in the 
complex A:-plane. The center-to-center seperation is i? = 6a and the wave number is measured 
in units of the inverse sphere radius a. 





0 5 10 15 20 

Re(k) [1/a] 


Figure 3: Quantum-mechanically calculated Eig- and Si^-resonances of the two-sphere sys¬ 
tem in the complex A:-plane. 
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Figure 4: Quantum-mechanically calculated E 2 g- and -F 2 n-i'esonances of the two-sphere 
system in the complex /c-plane. 





Figure 5; The three-sphere-configuration consisting of three equal spheres. The global and 
the local coordinate system of the first sphere are also shown. The remaining two local 
coordinate systems are obtained by rotations through ^ (sphere two) and ^ (sphere three) 
around the origin of the global system. All z-axes are perpendicular to the figure plane. 



Figure 6: Quantum-mechanically calculated Ai- and Ai^'-resonances of the three-sphere- 
system in the complex /c-plane. 







Figure 7; The four-sphere-configuration consists of four equal spheres at the corners of a 
regular tetrahedron (thick lines). In order to have a better visualization of the symmetries 
the tetrahedron is placed in a cube in whose center lies the origin of the global coordinate 
system {x,y, z). 
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Figure 8: Quantum-mechanically calculated ^i-resonances of the four-sphere system in the 
complex /c-plane. For numerical reasons only the first subleading resonances have been cal¬ 
culated. 



Figure 9: Comparison of the quantum-mechanically calculated resonances in the complex 
/c-plane of the totally symmetric representations of the two-sphere and two-disk systems. 
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Figure 10: Comparison of the quantum-mechanically calculated resonances in the complex 
/c-plane of the totally symmetric representations of the three-sphere and three-disk systems. 
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Figure 11: Comparison of quantum-mechanically (q.m.) and semiclassically (s.cl.) calcu¬ 
lated resonances of the completely symmetric representation of the t'wo-sphere system in the 
complex /c-plane. 






